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ABSTRACT: We consider the anomalous dimension of the fermionic highest states Tr tp L in 
the su(l|l) sector of M = 4 SYM at strong coupling. In the thermodynamical L — » oo limit 
it is described by a BES-like integral equation recently proposed by Rej, Staudacher and 
Zieme. The strong coupling regime of this equation is analyzed numerically and analyti- 
cally by Neumann expansion methods which have been developed for the s((2) sector. We 
compute analytically the first two terms of the strong coupling expansion and present nu- 
merical results for the next correction. We illustrate various specific features valid for the 
5u(l| 1) sector. In particular, at next-to-leading order, we find and solve a singular integral 
equation describing the scaling continuum limit of the Neumann coefficients. 
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1. Introduction 

The planar maximally supersymmetric JV = 4 SYM theory has an intriguing integrability 
governing the renormalization mixing of its composite operators [Q, ^ ^p. The associated 
integrable Hamiltonian is the dilatation generator of psu(2, 2|4) whose eigenvalues are the 
scaling dimensions. By AdS/ CFT duality [§, [|, its integrability properties are related 
to that of those of the dual type IIB superstring on AdS$ x S 5 [^]. 

All order Bethe Ansatz equations connecting the two sides of the correspondence 
have been proposed in [|8p. It is well known that these equations are asymptotic. The 
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weak coupling perturbative expansion of anomalous dimensions is correct up to a certain 
wrapping order depending on the length L of the operator that one is considering. In 
the thermodynamical limit L — > oo it is reasonable to assume that wrapping does not 
occur and it is possible to derive integral equations || [14|] describing exactly the flow of 
anomalous dimensions from weak to strong gauge coupling g defined as 

16tt z 

where A is the 't Hooft planar coupling. 

The Bethe equations describe integrable elastic scattering of internal elementary exci- 
tations in terms of 2 — > 2 processes governed by a specific S-matrix [jig]. The S-matrix is 
partially determined by psu(2, 2|4) symmetry Jl6| ] which completely fixes it up to a scalar 
factor, the dressing phase. It must obey the crossing relation of [|l3|]. The strong coupling 
expansion of the dressing phase can be matched to string theory calculations leading to 
the AFS leading approximation [10], and to the next-to-leading correction [O]. The crucial 



results of [ |T2| ] allowed to obtain the current complete form of the dressing phase proposed 
by Beisert, Eden, and Staudacher (BES) [p. 

The BES dressing factor is a clever analytical continuation of the known information 
on the string side of AdS/ CFT down to the weak coupling regime. If correct, it permits 
in principle to extract the perturbative weak-strong coupling expansions on the two sides 
from the analysis of the Bethe Ansatz equations bypassing the painful high-order stan- 
dard diagrammatic expansions. 

Recently, the BES formulation has been hardly tested looking at twist operators in the 
sl(2) subsector of N = 4 SYM. For these operators, one can consider the large spin N 
behavior of the quantum anomalous dimension |l8| ] 

Aat(<7) = /(<?) log N + O(N ). (1.2) 

The so-called scaling function f(g) is related to the cusp anomalous dimension introduced 
in l ^g , ^C| ] and related to the universal soft gauge boson emission. It is known up to three 
loops in the gauge theory [p]]. At four loops, it is known numerically with high preci- 
sion [^2|]. These results are in full agreement with the BES Bethe Ansatz calculation 
Remarkably, the weak coupling expansion of the BES integral equation can be extended 
to high orders with minor effort. 

At strong coupling, the expansion of the scaling function can be performed on the 
string side and the calculations of p3| ] leads to the expression 



«.)=».-^-AV, (1.3) 

7T 2 7T Z g 



where K is Catalan's constant. We write the explicit form of this strong coupling expan- 
sion to emphasize that it appears to be an expansion in integer inverse powers of g. 

Surprisingly, the analytic derivation of this expansion is rather non-trivial for reasons 
explained in pi}]. A direct numerical analysis of the BES equation has been proposed 
in [ p5| ] and led to a clear numerical approximation of the first three coefficient of the strong 
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coupling expansion in agreement with Eq. ( |1.3[ ). The analytical derivation of the first 
coefficient has been performed in [ p6| , |27| ] starting from the BES equation. The analysis 
is greatly simplified if the BES dressing phase is replaced by its leading strong coupling 
term, the AFS phase p8|]. Honestly this is not totally satisfactory since the main concern 
about the BES phase is its correct analytic continuation from weak to strong coupling and 
one would like to recover the AFS phase back from the BES phase. A solution to this 
technical problem has been proposed in the paper [24] which also deals with the other 
rank-1 su(2) and su(l|l) closed subsectors. 

If the aim is that of deriving the expansion of Eq. (fO| ) by efficient integrability tech- 
niques, then these investigations definitely suggested to analyze the strong coupling ex- 
pansion using the string Bethe equations leaving aside the BES equation. This approach 



was successfully adopted in [29] where the next-to-leading coefficient in Eq. ( [L3| ) was 
computed. The same result has also been obtained in [ f30| ] in the framework of the Baxter 
equation. 

However, quite surprisingly, this intriguing story is now basically concluded by going 
back to the BES formulation. In the beautiful paper [31], Basso and Korchemsky showed 
that the full strong coupling expansion of the scaling function can be recovered by an 
improved analysis of the BES equation. The new ingredient are crucial scaling assump- 
tions inspired by the numerical solution obtained by the methods of Recently these 
assumptions have been put on rigorous footings in [|^]. 

Given this rather rich landscape of methods and results, we believe that it would be 
interesting to adopt the same viewpoint in the case of other sectors of M = 4 SYM. Here, 
we consider the su(l|l) sector (see [|9p for its near pp-wave spectrum) and its maximally 
excited states of the form 

L = Trip L , (1.4) 

where ip is a gaugino component. These operators are eigenstates of the dilatation op- 
erator for each length (odd) L. The specific problem of deriving the strong coupling ex- 
pansion of their anomalous dimensions Al has been treated in some details in [|33|]. One 
motivation of this study was to understand if the Gubser-Klebanov-Polyakov A ~ A 1 / 4 
law (GKP) [|34|] is correct for these operators. A numerical and analytical analysis of the 
string Bethe equations with AFS dressing p7| ] (confirmed by the later calculation of [^4|]) 
proved that the GKP law is correct and led to the prediction 

lim ^ = 4s A 1/4 + • • • = + • • • • (1.5) 

The same result has also been obtained in [|^] using the light-cone Bethe equations pro- 



posed in [ p5[ ] in the general case and in [ p6[ ] for the su(l| 1) sector. 

These methods could not go beyond the leading term for various technical reasons. 
The only result at next-to-leading order is reported in [p6|]. An analysis of the quantum 
string Bethe Ansatz equations and of a su(l| 1) truncation of the AdS$ x S 5 superstring in 
the uniform light-cone gauge consistently gives the quantum anomalous dimension 
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Remarkably, the next-to-leading correction exactly cancels the classical contribution to the 
scaling dimension. 

Hence, an interesting question is that of recovering and possibly improve the above 
prediction starting from the BES integral equation describing this sector and derived in 
[|l4p. A first approach in this direction is the work p^ ] which however does not include 
the dressing effects in the strong coupling analysis. In this paper, we mimic the sl(2) 
analysis and test the direct methods of [f2|] and [^] to see if they are able to provide 
further information of the next-to-leading terms as well as the functional form of the Bethe 
roots distribution. In particular, we want to understand which is the origin of the A ~ A 1 / 4 
scaling compared to the A ~ \f\ regime which characterizes the 5 1(2) sector. We believe 



that this first step is mandatory if we want to apply the powerful all-order methods of [ pl| ] 



and [ p2[ ] to this sector. 

The detailed plan of the paper is the following. In Sec. (§) we recall the BES equation 
for the highest states in the su(l|l) sector in the thermodynamical limit. In Sec. @ we 
derive a linear problem for the coefficient of the Neumann expansion of the solution of 
the BES equation. In Sec. (||) we present a NLO approximation of the exact linear prob- 
lem. Sec. (|5]) is devoted to our numerical study of the exact and NLO linear problems. 
Finally, Sec. @ is devoted to the analytical proof of the NLO expansion of the anomalous 
dimension at strong coupling, based on the numerical insight described in Sec. (||). A few 
appendices are devoted to various technical problems and proofs. 

2. The BES integral equation for the su(l|l) sector 

The integral equation for this sector has been derived in [|||] and reads 

p[t) = e~* | J (2 g t) - 2 g 2 t J™ dt 1 [K m (2 g t, 2 g t') + 4 K A {2 g t, 2 g t')] p(i') J , (2.1) 

where p(t) is essentially the Fourier transform of the first-level Bethe roots. The explicit 
expression of the main kernel K m is 

h{t) j (O- Mt) Ji(f) 



K m (t,t') = K (t,t') + K x [tJ) 



t-V 



J*-o(t,t ) = 2 _ , 2 , 

The various terms admit the following Neumann expansion in series of Bessel functions 
*o(f,0 = ^ 7 ^(2n-l)J 2n _ 1 (t)J 2n _i(^) J 

n>l 

Kx{t,t') = ^2nJ 2n (t) J 2n {t'), (2.3) 

n>l 
n>l 
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The dressing kernel is 



/ — 'e*"-l 



K d (t,t') = 4:g 2 I dt"X 1 (t 1 2 5 i")^-r^o(25i",0. ( 2 - 4 ) 



and admits the Neumann expansion 

MM') = E E 8n ^7 " ^ ^n,2m-l J 2 »ft) hm-l{t'), (2.5) 
n>l m>l 

where the coupling dependent matrix Z n ^ m (g) is 

7 ffA f°°dt J n {2gt)J m {2gt) 

Zn,m(9) = J o J -TTi • (2-6) 

Unfortunately, as far as we know, this integral cannot be computed in closed form. Its 
accurate numerical evaluation is discussed in App. (|A|). 

Given the solution of the integral equation Eq. fl2.ip , the anomalous dimension of 
the su(l|l) operator O = Jitp L divided by the operator length can be computed in the 
thermodynamical limit as 

lim y= 2g 2 E(g), (2.7) 

L— too L 

where 

E(g) = A rdt J -^>p{t). (2.8) 
Jo 2gt 

The full anomalous dimension is obtained by adding the classical +3/2 contributions, the 
bare dimension of Tr ip L divided by L. 

2.1 Weak coupling expansion of the integral equation 

The perturbative expansion of p(t) for small g can easily be obtained as discussed in [fl4|]. 
The initial terms are 

p{t) = e-*{l -i(l + i) 5 2 + ^t{t 3 + 2t 2 -2* + 16) 5 4 + (2.9) 

-— t (t 5 + 3 £ 4 - 6 1 3 + 78 1 2 - 180 1 + 288 & t + 1044) g 6 + ■■■}, 
36 

and the associated energy E(g) reads 



E{g) = -4(! + = 2 " 8 9 2 + 58 5 4 - (518 + 32 £ 3 ) g 6 + • • • . (2.10) 

9 



3. The equivalent linear problem 

We introduce the new function a(t) = e* p(t) obeying the equation 

roo 

a{t) = J {2gt)-2g 2 t / d£ [K m {2gt,2gt') + AK d {2gt,2gt')] e"* a(t'), (3.1) 

Jo 
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and expand it in Neumann series 

a{t) = Y J Jn{2gt)a n) a = 1. (3.2) 

n>0 

The anomalous dimension is simply 

lim - = -4 5 <T!. (3.3) 

L—^oo Jj 

The integral equation can be written 

Y, W g t) a n = J (2 g t) - £ ^ m J m (2 5 1) W m , n a n (3.4) 

n>0 m>l n>0 

~ 16 X] Yl m ( 2k ~ 1 ) J 2m(2gt)Z2m,2k-lW 2 k-l,n^n, 

where 



W n>m {g) = / - e"*/^) J n (t) J m (t). (3.5) 
This integral can be evaluated in closed form with the result 

W (n)-a n + m r(n + m) x (3 6) 

w n , m {g)-g r(n + 1)r(m + 1) x V- b > 

x4i?3 ( \{n + m),\{n + m + l),\{n + m + l),\{n + m + 2) ._ 16 2 \ 
I n + l,m + l,n + m + 1 ' y 

We can split Eq. ( |3.4[ ) into the following two equations valid for m > 1 

0-2 m = -2m ^ V7 2m ,nO-n - 16m ^(2fc - l)Z 2m ,2k-i W 2 k-i,nCr n , (3.7) 

n>0 fc>l 
n>0 

0"2m-l = ~(2m- 1) ^W 2m -l,„Cr„. 

n>0 

The second equation can be replaced in the first obtaining our final form for the linear 
system equivalent to the BES integral equation 

Exact linear problem 



02m = -2m Y W 2m,nCT n + 16m ^2 2m ,2fe-lff2H, (3.8) 

n>0 fc>l 

C"2m-1 = ~(2m -1) ^^m-l.nffn. 



n>0 



3.1 Weak coupling expansion of the equivalent linear problem 



As a check we can recover from the linear problem Eq. ( p.8[ ) the weak coupling expansion 
of the energy. We start from the easy expansions 

Z n , m = 0(g n+m ), W n<m = 0(g n+m ), (3.9) 
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finding for instance 

Z 2l i = CsS 3 -10C55 5 + 105Cr5 7 + ••• , Wi,o = 5-3 5 3 + 20 5 5 + -- - , 

Z 2 ,s = 2t 5 g 5 -35t 7 g 7 + ■■■ , W 1A = g 2 - 6 5 4 + 50 g e + • • • , (3.10) 

Z 2 , 5 = 3Cr<7 7 + -- - , Wi,2 = 3 3 -10 5 5 + 105 5 7 + -- - . 

Then, we expand the linear system with the position 

p>0 p>0 

under the condition <7o = 1. We immediately find 
<ri = -5 + 4 5 3 -29 5 5 + -- - , 

<r 2 = -g 2 + (10- 16 Cs)<? 4 + (-100 + 64(3 + 160 c 5 ) g 6 + ••• , (3-12) 

as = - 5 3 + 18<? 5 - 240 5 7 + -- - , 

a 4 = -<7 4 + (28 - 32 C 5 )5 6 + (-482 + 128C5 + 672C 7 )s 8 + --- , 



which reconstructs the previous expression of a(t) order by order in g 2 . 
3.2 Basso-Korchemsky zero-mode formulation 



Let us derive zero-mode equations in the spirit of [31]. We split 



a{2gt) = a + {2gt) + a_{2gt), 
a + {2gt) = J {2gt)+J2 J 2n{2gt)a 2n , (3.13) 



n>l 



a_{2gt) = J2 J 2n-i(2gt) 



n>l 



If n and m are positive and have the same parity then we have 

r°° dt 1 

/ - Jn{t) Jmit) = — 8 n , m . (3.14) 
JO t 2,71 

From this relation and the definitions of W and Z it is straightforward to derive the fol- 
lowing equations (n > in the first one) 

l ~hn{r) [a + (r) + -e^\a + {r) + g _(r)) - _ i a_(r)| = 0, 



y J 2 n+i(r) {a_ + V^Wr) + a_(r)) r = 0. (3.15) 



Hence, 



a+(T) + ^ e -^ 2 ^(a + (r) + a_( T )) - _ l *_(t) = x+(r), (3.16) 

a _ + I B -r/(2 fl ) (ff+(r) +CT (r)) = x _ (r)) 



-7- 



where X±( r ) are the most general solutions of the equations 

f°° dr 

/ — hn{r) X+(r) = 0, (n>0) (3.17) 

JO T 

f°° dr 

/ — J 2 n + i(r) X-(r) = 0. (3.18) 

As in the case of the sl(2) sector, these zero-mode equations admit many solutions and 
one must find a way to fix their contributions. In this paper, we focus on the next-to- 
leading solution of the integral equation and shall not discuss in details the zero-mode 
equations. Notice that this is already a non-trivial task since the non-homogeneous piece 
in the integral equation gives only c+(0) = 1 and it can be seen that the normalization of 
£7_ is not fixed by the zero-mode equations. It must be found from the integral equation, 
or the equivalent linear problem. 



4. Next-to-leading order linear problem 



Following the approach in it will be very convenient to study an approximation 
to the exact linear problem Eq. ( f3!8[ ). This is obtained by replacing the integrals Z and 
W by the first terms of their asymptotic expansions. We choose the following level of 
approximation 



where 



and 



Z n , m (g) = gzZ >m --wZ tm + 
W„, m (5) = < m + ---, 



Zj 2m,2k-l 



1 

4m 



2m - 1' 
1 

2m + 1' 



m 



k = m + 1 



< m (s) 



dt 



Jn(t)Jm(t) 



1 

2n' 

2 sin((n - m)f ) 



7T 



m z 



n 



m 



(4.1) 



(4.2) 



(4.3) 



This approximation will be denoted as next-to-leading order (NLO). The reason is that 
the strong coupling expansion of Z and W goes in inverse powers of g times logarithms, 
whereas the natural expansion parameter for the anomalous dimension is supposed to be 
l/y/g, with possible logarithmic corrections. 

Thus, at NLO, the exact linear problem Eq. ( ^ ) can be written in the following form 
(m > 1 in the first equation and m > in the second one) 
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NLO linear problem 



- a 2 m = - 10 m £ W2m,2 n+l °2 n+l + 4 Q [ j^—j + ) , (4.4) 



n>0 



2°2m+l = -{2m + I) J2 W 2m+l,2n^2n- 



n>0 



Before attempting an analytical study of this problem, let us present a numerical in- 
vestigation of the exact and NLO linear problems. 

5. Numerical analysis of the linear problems 

5.1 Strong coupling expansion of the anomalous dimensions 

As a first step, we have studied numerically the exact linear problem Eq. ( ^) . To this aim, 



we have truncated the mode expansion up to N modes and have studied the convergence 
of the solution with increasing N at various couplings g. Examples are shown in Fig. (||). 

A better behavior is obtained by studying in the same way the NLO linear problem 
Eq. (fO|). This is nice as illustrated in Fig. (||), but what is the accuracy of the NLO approx- 
imation ? Reasonably, it should work well for large g at least with the aim of reproducing 
the NLO expression of the strong coupling expansion. This assumption is explicitly tested 
in Fig. (^). As the figure shows, the agreement between the exact and NLO calculations is 
quite good in the considered range. 

Given these positive results, we have attempted to extract the strong coupling expan- 
sion of the energy density by a numerical fit of the NLO solution at quite large values of 
g. The analysis is shown in Fig. (B). The leading term ~ g 1 / 2 is evident. The functional 
form that we have chosen and which is rigorously suitable for the NLO problem is 

A , 3 log g 1 

lim — = V27Tgci - -c 2 H —c 3 + — C4-\ . (5.1) 

L^oo L 2 ^Jg ^Jg 

The best fit gives 

\d - 1|, |c 2 - 1| ~ 10- 5 . (5.2) 

strongly supporting the prediction in p6| ] 

A3 
lim — = ^/2¥g~ -- + ■■■ . (5.3) 

L— >oo h A 

In particular, we have recovered the intriguing cancellation of the classical +3/2 contribu- 
tion. 

The next corrections are definitely beyond the reach of the NLO approximation and 



must be extracted from the solution of the exact linear problem Eq. ( |3.8[ ), although the ^fg 
dependence requires to push the exact calculation to quite large g. We have worked up to 
g = 180 which is rather challenging from the computational point of view. The solution 
of the exact linear problem is well fitted by the following terms 

,. A , 3 0.12(1) log + 0.28(1) , x 

lim- = ^2^-- + K - L -^j=- + (5.4) 

L->oo L 2 *Jg 
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as illustrated in Fig. (g). 

The agreement is rather remarkable, although we are using a simple two-parameter 
fit. Of course, the logarithmic NNLO term deserves special attention. For instance, it ap- 
pears to be absent in the sl(2) sector. Therefore, despite the good numerical evidence, we 
prefer to be conservative and we do not exclude the possibility that it could be a numerical 
fake. For instance, it could be due to fitting over a too narrow range of values of the cou- 
pling g. As another possibility, it could mimick a non-perturbative contribution once the 
log-free perturbative tail is subtracted. In conclusion, we believe that additional numeri- 
cal study at larger g would be important to assess the NNLO corrections in Eq. ( fT4| ) in a 
firm and safe way. This is a worthwhile issue in order to exclude possible order-of-limits 
problems in comparing semiclassical string theory calculations against BES-like integral 
equations. 

5.2 Shape analysis 

It is interesting to analyze the shape of the solution of the NLO problem at various g. This 
can be done by plotting at increasing N the even and odd parts of a n after multiplication 
by a trivial alternating sign, i.e. 

even part : (-l) n cr 2 „, (5.5) 
odd part : (-l) n a 2n +i- 

The result is shown in Figs. (£>[ g, g) for g = 10, 100, 1000. 

The conclusion is that the even/ odd part of a n tends to smooth curves (apart from 
the alternating sign) and a non-trivial scaling shape emerges. This is an important issue 
that we can partially analyze quantitatively. 

Indeed, from the numerics, it is natural to set 

A a W 

cT 2m+ i = (-l) m (2m+l)- + ^ + ... , (5.6) 

\/9 9 

a 2m = a^ + 0(g-^). 
Replacing this Ansatz into the second equation of Eq. ([O]) and using = 1 we find 

(2 m + 1)2 + ^(2m + l) 2 -(2n) 2 2n K ' 

A simple solution, again well supported by the numerics, is 

afi = 2(-ir. (5.8) 

The normalization of this contribution is fixed and this is not surprising since it is related 
to the normalization cr + (0) = 1 in the zero-mode formulation. Instead, things are quite 
different for the odd part. 
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Replacing Eq. ( |5~8| ) into the first equation of Eq. (fOJ) we get 



-•2(-l) m =4 ^ ±L _^2rn- L (g 9) 

2 v 7 2m+ 1 2m - 1 v 7 



The general solution is 



4m+i = (-!) m (2m+ 1) f|(2m+ 1) + a ) , a, constant. (5.10) 



It depends on an arbitrary constant a which is not fixed in this naive expansion. In con- 
clusion, we have obtained the following expansion 



a 2m = 2{-l) m + ■■■ , (5.11) 
<J2m+i = (-l) m (2m+l) 



A + Ig (2m + 1) + a ) + ... 

where A and a are not determined. 

This indetermination is clearly related to the zero-mode contributions in the zero- 
mode formulation. Summing at this order the functions ct±{t) (r = 2gt) we find the 
results 

ct + (t) = cos r H , (5.12) 

1 T , ^ ( A a\ 3 
a_(T) = -TJ (r) (^ + ^+T^ C0S - + ---. 

which satisfy the leading order zero mode equations. The first term in the second equa- 
tion contains the zero-mode %_ = r Jo(t) with a basically undetermined coefficient plus 
a contribution which is inherited from the even leading zero mode x+ = cos r whose 
normalization is fixed by cr + (0) = 1. 



5.3 Scaling 

We now ask ourselves the basic question whether it is possible to derive analytically the 
results 




(5.13) 



which are required to reproduce the observed strong coupling expansion of the anoma- 
lous dimension at next-to-leading order. 

A simple way to achieve this result is based on the above remarks about scaling. Our 
expansion for the odd subsequence 02 m+i can be written in the following suggestive form 

f ir<T _ 9 ,™ + 5 3 /m + i\ 2 2a m+| 
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This means that if we take the limit g — > oo with fixed x = —jf- we find 

(-l) m V2m+l = F{X) = F$\x) + 4= F^\x) + ■ ■ ■ > 

F^\x) = 2Ax + ^x 2 + ••• , 
F^\x) = 2ax + --- . 



(5.15) 



This result is numerically tested in Fig. @. Hence, the desired coefficients A and a can be 
extracted as 



A = - lim v ' 



1 ,. FM(x) 
a = - lim 



(5.16) 



2 x^o 2x 2 x^o 2i 

We now illustrate a way to compute the functions and F^ 1 ' with the required accu- 
racy. 



6. The shape equation in the continuum limit 

Let us define for n > 

<7 2n = (-l)"z+, a 2n+1 = (-l)"z n . 
The NLO linear problem can be written 
2 1 



3 + 
2 Zm 



^ 7r(2n + l) 2 -(2m) 2 n 5 \2m+ 1 2m- 1 



^z m = (2m + l) £ 2 



n>0 



7r(2n) 2 - {2m + 1) 



z+, m > 0. 



2 n 



The second equation can be inverted by using the result 



E 



7T 



* m>0 + ^ ~ ( 2?? ) 2 ( 2m + X ) 2 ~ ( 2n ) 2 



16 p 2 
32p 



, n = 0, 
2 , n > 0, 



(6.1) 



m > 1, 



(6.2) 



(6.3) 



The exchange of summation is safe due to the asymptotic properties of z n as can be 
checked a posteriori on the numerical solution. We obtain 



P * ^ ( 2 P) 2 -( 2m + 1 ) 2 2m + 1' 



(6.4) 



Replacing in the first equation of Eq. ( |6.2| ) and setting 



2n + 1 ! 



(6.5) 



-12- 



we arrive at the single equation 

8 \_^_^T t 40 a V" 2m + 1 & 

3 " " 7r^ (2^-(2m + l)2 U+ Stt ^ (2^-(2m + l)2 U 

(6-6) 

where a = 1. The reason for introducing this parameter is that it will be irrelevant to the 
aim of computing the constant A. Therefore we shall be able to choose it in an optimal 
way. At the numerical level, we can solve Eq. ( |6lj| ) by imposing as usual the condition 
£w = for some large N to be increased up to convergence. In Fig. ( |T0| ) we show the 
solution for four values of a (including the original a = 1) at g = 100 with N = 600 
modes. As one can see, the curve for £ n starts at the same value. 

By the way, we immediately see that the determination of the constant A in Eq. Q5.13 ) 



is non trivial. The problem Eq. fl6.6D is solved with the boundary condition £ n — i> for 
n — » oo. This position produces a unique solution and a non ambiguous prediction for £o 
which depends non-locally on the whole solution. 

The rigorous proof of independence on a is illustrated in Appendix @ where it is 
shown to hold at least in the range < a < 9/5 which includes the default value a = 1. 
Due to this remarkable feature we can make the easiest choice a = and consider the 
simplified problem 

O AO _2 

3»({.-iH)- 7 £ w .; 2 „ + 1) , U-2 = .. (6.7) 

We repeat that this problem is well posed in the following sense: We truncate the sequence 
{£m} to < m < N and impose £w = 0. The above linear problem has a unique solution 
with a definite point- wise limit as N — > oo. 

6.1 Continuum limit 

Based on the numerical analysis, we define the continuum limit of Eq. (|6.7| ) by setting 
x = vl \[9 an d assuming the following scaling form valid at fixed x and large g 



P + l/ 2 

V9 J V 



Cp = — f\^7=-)+---, (6-8) 



where / is a suitable smooth function. The half-integer shift in the argument of / is im- 
portant. Indeed, the results of App. Q imply that 

(i) f(x) obeys the following Cauchy integro-differential shape equation 

9 r°° x 2 3 
fV-Tti x^ f M dy -l = > (6.9) 
/(+oo) = 0, 

The relation between f(x) and F^°\x) defined in Eq. ( |5.15[ ) is 



fix) = — , (6.10) 
2 x 
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and therefore A = £(0). Notice also that the shape equation immediately implies 
£'(0) = | under very mild assumptions on its solution. This is in agreement with 
the second term in F^°\x). The solution of the shape equation is invariant under 
constant shifts £ — s> £ + constant. This freedom is fixed by the condition £(+oo) = 0. 
An equivalent form of the shape equation is easily obtained by manipulating the 
principal integral arriving at 

9 r°° v 2 

f'W-^zt ^72/(2/)^ = °' 

2 7T Jo £ — y 

(6.11) 

roo 



r°° 7r 

J f{x)dx = --, /(+oo) = 0, 



where now the (homogeneous) first equation has a rescaling freedom which is fixed 
by the area constraint in the second line. 



(ii) the corrections to Eq. ( |678[ ) are 0(1/ g), i.e. 

1 { <n 4- \ l0.\ 

(6.12) 



+ 



This immediately implies that F^\x) = and therefore a = which is the second 
part of Eq. ( gig ). 

To summarize, we have shown that the NLO expansion of the anomalous dimension 
can be proved analytically if we are able to show that the solution of the shape equation 
Eq. |x9|) satisfies 

/(0) = -\J\ ■ (6-13) 

This will be proved in the next section. Later, we shall derive additional more complete 
information about the function f(x). 

6.2 The value /(0) = A 

The value /(0) can be computed even without knowing the complete solution. To this 
aim, we write 



>'« = h f 3^ * = I f 2 » 2 + ^;/ 2 >;^-» 2 > ,(,) d y 
2 7T Jo x — y 4 7rJo £ — y 



roo 



9 f 00 4- v 2 9 

" ' f^f(y) d y-— I f(x)dx. (6.14) 



47T Jo x 2 — y 2 47T 7o 

Using the known area constraint in Eq. fl6.1ip , we find 



9 f°° x 2 + y 2 3 



/'(«) = 7T f 7^2/(y)^ + «- ( 6 - 15 ) 

4 7T Jo 2: 2/ 8 
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Hence, using /(+oo) = 0, we get 

^oo ^oo 

/(0) 2 = - / (/(a) 2 )' dx = -2 f{x) f'(x) dx = 
Jo Jo 

Q poo poo ™2 I y2 a poo 

= -— / dxl d y—^f(x)f(y)-- f(x)dx. (6.16) 
2 7T Jo Jo % — y 4 Jo 

The first term vanished by antisymmetry x <-> y. The exchange of a standard integral and 
a principal value integral is allowed by the results of [[40|]. Using again the area constraint, 
we find 



4 V 67 8' 

As a comment, we remark that this result can also be obtained directly from the discrete 
equation Eq. (|6^). The method is analogous to the above derivation. We multiply by 
£ p + £ p -i and take a sum over p. After symmetrization of the double sum, one easily 
proves the above boundary value as the truncation number N — > oo. 

6.3 Complete solution of the shape equation 

A canonical form is obtained by scaling z = x and setting / zj = h(z). It reads 

1 f°° w 2 
h'(z)--f — =h{w)dw = 0, 

7T Jo Z 2 - W 2 

(6.18) 

roo ^ 

^ h(z)dz = -^=, h(+OO) = 0. 

To solve the shape equation in canonical form we follow the methods of p3|]. We change 
variable and define h^-Ju) = g(u). In the new variable u, the equation reduces to 

1 r°° \ft 

g'(u)-—i -?—g{ut)dt. (6.19) 

4 7T Jo 1 - t 

The Mellin transform [M] 

poo 

g(s) = / u s ^g(u)du, (6.20) 
Jo 

has fundamental strip (0,1) and obeys the difference equation 

m 



9(s ~ 1) 



4(s- 1) cot its. (6.21) 



The solution which is associated with a g(u) regular in u = and vanishing for u — > oo is 
easily shown to be 

g(s) = 4 s T(s) $(-s) /c, (6.22) 

where k is a constant to be determined by the area constraint (or equivalently by the 
condition on /'(0)). The function §(s) is the ratio 

, x G(l - s)G(| + s) 
*( fl ) = — r> 6-23 



-15- 



where G(z) is the Barnes G-function discussed for instance in [42] and defined as 

W) I I ( 1 + - f e- z+ f 

2 K '2 



G(z + l) = (27r) z / 2 exp(-lz(z + l)-J 7 EZ 2 ) Y\(l + ^) e 



n>l 



It obeys in particular 



G(z + 1) = r(z)G(z), 
G(l) = 1, 



(6.24) 



(6.25) 



from which it is easy to prove that the difference equation is satisfied. Useful properties 
of the ratio <t>(s) are the infinite product expression 



$( a ) = v^f n 



1 + 



the functional relations 



(6.26) 



$(s) = <!>(— s — -) = "^(s — 1) tan 7rs, 



and the expansion in 5 = 



1?7r 2\ .3 



3T + - 



(6.27) 



(6.28) 



From these relations we can look at the poles for negative half-integer 5 and for posi- 
tive integer s and use the known properties of the Mellin transform to derive asymptotic 
expansions of g(u) around u = and u = +00. The result is 



k~ 1 g(u) = A - u 1 ' 2 + v?l 2 ( -±- log u + 11 37E 



67T 



18 



+ G(u 5/2 log 2 u), (6.29) 



« x 5(u) 



4 1 16 1 

+ 



l-KU 7T 3 / 2 U 2 

Taking k = -y/2/4 we recover the correct expansion 



(- log u + 2 log 2 - 7s ) + 0(u~ 3 log 2 u). 



/(*) 



/F 3 

+ -5C + 

4 



(6.30) 



As an additional check of the solution, we can evaluate the area constraint finding again 
agreement 



jf M Z )^ = ijQ)=-^,Hm4-r W *(- s ) 



(6.31) 



f .2.^/2)^=-^. 



A numerical comparison of the above asymptotic expansions with the numerical solution 
of the shape equation is illustrated in Fig. (11). 
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7. Conclusions 



In summary, we have considered the BES-like equation describing the anomalous dimen- 
sion of highest states Tr ip L in the su(l|l) sector for L — > oo proposed in [14]. We have 



followed the approach pioneered in pq , [26| ] for the st(2) twist operators. It is based on 
the Neumann expansion of the Fourier transformed Bethe root density. Although many 
technical differences arise, we have been able to analytically derive the next-to-leading 
order expansion of the anomalous dimension at strong coupling. Some interesting infor- 
mation on the NNLO terms is also presented based on numerical calculations only. Our 
results reproduce at leading order the result already discussed in the literature, although 
with different methods. At next-to-leading order, we find agreement with the proposal in 



[p6|]. Besides, several specific interesting features appear. In particular, at NLO we have 
found and solved a singular integral equation describing the scaling continuum limit of 
the Neumann coefficients. 
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A. Numerical evaluation of Z n>m (g) 

A simple strategy to numerically evaluate Z n ^ m (g) with arbitrary precision is based on the 
exact splitting 

N r°° H+ r°° Hf P~(. N + 1 ) t /( 2 9) 

Z n , m (9)=T,J J^ n+1)tf{29) Ut)Jm{t) + J J ( , t/(2g) _ 1 Ut)J m (t). (Al) 

71 = ^ ^ 

The first finite sum can be computed in terms of the known function W U)rn {g) since 

r *-(«+!)*/(**> Ut) Jm( t) = Wn , m . (A.2) 

Jo t \n + lj 

The second integral can be computed by splitting the integration interval in subintervals 
separated by the roots of the product J n (t) J m {t)- The resulting sub-integrals build an 
alternating sum whose terms can be computed numerically without any difficulty. The 
optimal value of N is N ~ g, but of course the result is independent on N. 

B. On the a-independence 

The shape equation for a = is Eq. Q6.ll ) that we repeat for the reader's convenience 



9 f°° y 2 



2 7T Jo £C — y 
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If the parameter a is not zero, the same steps starting from Eq. 
with the simple substitution in the numerator of the integral 



V 2 -Pxy, 



6 = - a. 
y 9 



give a similar result 



(B.2) 



We now show that for < j3 < 1 the first two terms in the asymptotic expansion of f(x) 
around x = do not depend on a. 

To this aim and with the notation of Sec. (6^3) we immediately find that Eq. (|6.21 ) is 
replaced by 

5 < S > - -4( S -1) 



Setting 
we find 



g[s — 1) ~ v ~ tan its + @ cot its' 
g{s) = /c4 s r( S )3o(s), 
9o( s ) _ 1 



(B.3) 
(B.4) 
(B.5) 



ga(s — 1) tan tts +/3 cot 7rs' 

From now on, let us consider the case < @ < 1. The zeroes of the denominator of the 
right hand side are located at the positions 



n±ia, 



n e 



(B.6) 



a = — arctanhv^. 



7T 



Looking at the zeroes and poles (with their order) of a go which is regular in the funda- 
mental strip < s < 1 we find the solution 



3o(s) 



G(l-s) G(f-s) G(l + ia + s) G(l-ia + s) 
G(l + s) G(f + 5) G(l + ia-s) G{l-ia-s) 



(cosh 7ra) 



2s 



(B.7) 



where G{s) is the Barnes G-function introduced in Sec. ( J673| ). The analytic structure is 
manifest in the Weierstrass infinite product form which is 



3o(s) = v^cosh -waf 3 f\ 



n>l 



1 " 1 



"+| 1 + n+ia 



1 + 



1 + M + 

n 



T 1 



n+i a 



(B.8) 



The area constraint and the asymptotic expansion of g(s) have a a-dependence which is 
all contained in the factor 



G(l + za + S ) G(l-ta + 3) 25 

P(5,a) = -— ; r— ; r ( C0SI1 7T a ) . 

VK ' ' G(l + za- s) G(l-ia- s) K ' 
The area constraint on the total integral J °° f(x) dx is a-independent since 



(B.9) 



cosh ir a 



G(| + ia) G(| -ta) 
G(| + ta) G(| -ifl) 

r(^ + i a) r(- - i a) cosh ira = tt. 



(B.10) 
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The coefficients of the leading and subleading terms of f(x) for small x are also a-independent 
since they are proportional to 



p(0,a) = l, p^--,aj=-. (B.ll) 
Finally, the next correction does depend on a since 

p(-l,a) = -Ar tanh 2 (7ra) = -. (B.12) 
As a final remark, we notice that the a-independence of /(0) can be proved easily using 



the method of Sec. fl(x2|) , assuming always /(+oo) = 0. Indeed the a-term just add a xy 
contribution to the numerator in the double integral appearing in Eq. ( |6.16| ). Since this is 
symmetric under x <-> y, the proof valid for a = is completely unchanged. 

C. Detailed continuum limit of the shape equation 

Let us consider Eq. <^J§ with the position p = x <fg and £ p = - 4= / ( j + Jl 2 ^. We shall 
analyze the large g expansion at fixed x for the two terms of Eq. ( |6^ ) to next-to-leading 
accuracy. The finite difference is immediately expanded and gives 

9 & - = /'(*) + ^ A*) + 0( V5 2 )- (CI) 

The infinite sum is more complicated. We use the Euler-MacLaurin formula 

E = f a F{m)dm+ na)+ 2 m + E Jgj [l*"" 1 ^) " ^V)] , (C2) 



k> 



where 

It is convenient to split the infinite sum in order to easily display the appearance of a 
principal value at leading order. Thus, we write 

/ J3-1 OO \ 

S = Y F(m) = I E + E (C.4) 

m>0 \m=0 m=pj 

Applying the Euler-MacLaurin formula we get 

a = ( r 1 + n *» + f(o)+j?( v i)+j?w + (c.5) 



'0 */p 

where R p is the sum involving higher odd derivatives of F. The two integrals can be 
written 

+ jT) F W *, = (/^ + £\ ) f(y) * = (C.6) 

io r-f 4^/5 
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The other two terms in the right hand side of Eq. flC.g ) have the expansion 

F(0) + F(p - 1) + F(p) _ 2 /(0) + f{x) - 2 x f'(x) +0 (l\ (Q7) 



2 4^ V<7 

i? p = O ( ~ ] . (C.8) 



.9 

In conclusion, summing the various terms, we find 



S = £ F(m) = /°° /(y) + O (±) , (C.9) 



m>0 

with complete cancellation of the 1 / ^fg correction. 
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Figure 1: Convergence of the solution o\ of the exact linear problem Eq. (3.8) as the number of 
truncated modes N is increased. The convergence is clearly slowing down as the coupling in- 
creases. 
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Figure 2: Convergence of the solution of the NLO linear problem Eq. (4.4) as the number of trun- 
cated modes 7\f is increased. The convergence slowing down is not is real problem even at rather 
strong values of the coupling. 
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Figure 3: Comparison between the solutions of the exact linear problem Eq. (3.8) and its NLO 
approximation Eq. (4.4). For all values of the coupling, the number of truncated modes N has 
been increased up to convergence. 
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Figure 4: Analysis of the leading and subleading strong coupling expansion of the energy density 
from the numerical solution of the NLO linear problem Eq. ([O]). 
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Figure 5: Numerical fit of the NNLO correction to the strong coupling expansion of the anomalous 
dimension. The best fit values of the two parameters are a = 0.12(1) and ai = 0.28(1). These 
numbers are basically unchanged if the first data points (at smaller g) are removed from the fit. 
Also, the point with the largest g = 180 is very well reproduced by the best fit computed with the 
points with g < 180. This is a check of the proposed functional form. 
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Figure 6: Shape of the NLO solution at g = 10 as the mode cut N is increased. 
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Figure 7: Shape of the NLO solution at g = 100 as the mode cut N is increased. 
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Figure 8: Shape of the NLO solution at g = 1000 as the mode cut N is increased. 
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Figure 9: Scaling analysis. Odd part of the NLO solution at g = 1000. The numerical shape is 
compared with the analytical result Eq. ( [5.14 ) with the correct choice Eq. ( [5.13 ). 




-30- 




Figure 11: Comparison between the solution of the NLO linear problem at a = and the asymp- 
totic expansion of the shape equation at x = and +oo. The first curve is obtained with g = 100 
and N = 600. 
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